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Abstract. The asymmetry of D~ and D + meson production in 7r - iV scattering observed by the E791 
experiment is a typical phenomenon known as the leading particle effect in charm hadroproducton. We 
show that the phenomenon can be explained by the effect of light quark fragmentation into charmed hadrons 
(LQF). Meanwhile, the size of the LQF effect is estimated from data of the E791 experiment. A comparison 
is made with the estimate of the LQF effect from prompt like-sign dimuon rate in neutrino experiments. 
The influence of the LQF effect on the measurement of nucleon strange distribution asymmetry from 
charged current charm production processes is briefly discussed. 

PACS. 14.60.Pq - 12.15.Ff 



1 Introduction 

The leading particle effect in charm hadroproduction has 
been observed by many experiments 1 2 ,'! HE] - The main 
feature of this effect is the enhancement in the production 
of the charmed hadrons that carry the same valence fla- 
vor of the incident hadron in the forward region, i.e., at 
positive xf, where xf is the Feynman variable for the 
produced hadron, xf = p*/p* max ^ with P* being the 
momentum along the beam direction in the cm. frame 
of the colliding hadrons. The result of the E791 exper- 
iment [3] is a typical example of such phenomenon with 
high statistics. In the E791 experiment, a 500 GeV ir~(ud) 
beam is incident on a fixed target, and an obvious excess 
of D~ (cd) (having the same valence flavor d with 7r~) over 
D + (cd) has been observed. The asymmetry variable 
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da(D-) -da(D+) 
da{D-) + da(D+) 



vs. xp and P% is used to describe the effect, and a clear 
rise of the asymmetry A with xp has been observed. 

On the theoretical side, because perturbative QCD 
(PQCD) predicts no asymmetry at leading order (LO) 
and very small asymmetry at next to leading order (NLO) 
for charm and anti-charm quark productions 6,7, the ob- 
served asymmetry is generally attributed to the hadroniza- 
tion processes. The "beam drag effect" [8] implemented in 
the PYTH1A Monte Carlo explains the asymmetry through 



the color strings between the produced charmed quarks 
and the beam remnants that moving in the same gen- 
eral direction. The charmed hadron produced through the 
decay/collapse of the low mass string may have more en- 
ergy than the original charmed quark through the pull of 
the fast beam remnants. By adjusting some parameters 
of the model, they can reproduce the observed asymme- 
try. Another explanation is the intrinsic charm coalescence 
model [9 , in which an intrinsic charm quark of the projec- 
tile combines a valence quark of similar rapidity to form 
the charmed leading particle. However, this model predicts 
much smaller asymmetry than what is observed experi- 
mentally. A further work is the heavy quark recombination 
mechanism |10j . in which a produced heavy (anti)quark 
recombines a light parton of similar velocity that partic- 
ipates in the hard scattering process. The work employs 
a simple PQCD 0(al) picture and describes the high xf 
and Pj, distributions of the asymmetry very well. There 
are still other works of quark recombination and quark- 
gluon string model [TTTfT2] ■ 
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Although so much work has been done, we should still 
explore other possibilities. As first promoted by Dias de 
Deus and Duraes [13] in this field, we investigate the pos- 
sible contribution of light quark fragmentation (LQF) into 
charmed hadrons. The LQF effect is an old idea originally 
suggested by Godbole and Roy [14] to explain the unex- 
pected high rate of prompt like-sign dimuon production 
from many neutrino experiments [151116). Although there 
are deviations in experiments and later high statistic ex- 
periments tend to show a smaller effect than early exper- 
iments, the LQF effect could not be ruled out experimen- 
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tally. When the leading particle effect in charm hadropro- 
duction is observed, Dias de Deus and Duraes [13] point 
out that the generally neglected fragmentation, e.g. d — > 
D~(cd) could be a possible mechanism to explain the 
observed asymmetry. When the NuTeV anomaly [17] is 
shown to be hopefully settled by the nucleon strange asym- 
metry |18 P 19 P 20 P 21 P 22 P 23 P 24|. while on the other hand, 
CCFR and NuTeV measurements do not show evidence 
for the strangeness asymmetry [25ll26| . we point out that 



the size of the LQF effect is roughly proportional to the 
number of charmed sea above £, n c ~ c(x)dx. This will 
yield a E q dependence of the LQF effect, which is nearly 
a linear rise with E q . This trend is consistent with the ob- 
served prompt like-sign dimuon rates in experiments [15l 

ng. ' 

From the pick-up process, the outgoing light quark can 
hadronize into a charmed hadron with a fraction z of its 
initial energy along its initial direction in the nucleon rest 



the LQF effect is possible to influence such measurements [27] frame The flavor of the Hght quark ig mQst probable to 



In this paper, we will investigate the LQF effect in charm 
hadroproduction, and estimate the size of the LQF effect 
from the E791 experiment. In section II, we will describe 
the physical mechanism of the LQF effect, and present its 
parametrization. A detailed description of our calculation 
will be given in section III. The results will be shown and 
be discussed in section IV. In section V we will give the 
summery and conclusions. 



2 The LQF effect 

In our model of the LQF effect, we differentiate the favored 
fragmentation q — > D(cq) or q — > D(cq), which is non- 
negligible, from the unfavored fragmentation, e.g., q — > 
D(cq'), which could be an indirect effect from LQF due to 
charm flavor conservation and can be neglected, since it is 
mainly confined near target area, i.e., in large negative Xp, 
where experiments were intangible. In 7r _ iV scattering, for 
example, since there are more d quarks than d quarks in 
7r~, more d than d will participate in the hard scattering, 
and with the favored fragmentation of d — > D~(cd) and 
d — > D + (cd), more D~ than D + will be produced in the 
forward region, and this is in accord with what is observed 
in experiments. 

In this work, the picture of the LQF effect has some 
differences with works of other authors. We propose a 
pick-up mechanism of the light quark fragmentation. Af- 
ter hard scattering of subprocesses, the outgoing light 
quark will experience a period of hadronization. This pe- 
riod AT ~ 1/M is relatively longer compared to some 
processes with large energy transfer At ~ l/AE. Then, 
through large momentum transfer of the strong interac- 
tion, the outgoing light quark can pick up a charmed (anti- 
) quark from the nucleon sea, forming a cluster, which can 
then decay into a charmed hadron. 

Although in the pick-up process, the light quark gives 
most of its energy and momentum to the nucleon system, 
the Q 2 of the interaction can be very small along the initial 
quark direction, and thus the process is actually a non- 
perturbative one and is not largely suppressed by a s (Q 2 ). 
For a rough estimate, one can consider a massless outgo- 
ing light quark with energy E q in the nucleon rest frame. 
Through the pick-up process, the light quark turns to a 
constituent (with constituent mass to be about 0.3GeV) 
of the produced D meson along its initial direction. It 
is easy to get the squared momentum transfer Q 2 ~ 0.5 
GeV 2 . The lowest momentum fraction of charmed sea from 
the nucleon in this process should be £ ~ m c /AE q , and 



be a valence flavor of the produced hadron. 

For a quantitative estimate of the LQF effect, we con- 
struct the parametrization of the light quark fragmenta- 
tion function (z, E q ), which describes the fragmenta- 
tion of light quark d into D~ meson (including D~ from 
D*~ decay), with E q and z being the energy of the light 
quark and the energy fraction of the produced D ~ meson, 
z = Ejj/Eq, in the nucleon rest frame. We assume charge 
symmetry D°~(z,E q ) = (z,E q ) = D f {z,E q ) and 
take Df(z, E q ) — C(E q )- P(z), where P(z) is a normalized 
parametrization and we use the form P(z) oc z a (l~ z) [28] 
with a being a free parameter. The form of C(E q ) should 
be an approximate linear rise as we have mentioned. We 
take the form C(E q ) — a ■ (E q — Eq) when E q > Eq, where 
Eq corresponds to an energy threshold for the LQF effect. 
We also suppose a suppression for the fragmentation of 
low transverse momentum light quarks, which correspond 
to remote interactions. Thus we include the suppression 
factor (1 — m 2 H /p 2 T ) with the restriction p 2 T > m 2 Hl where 
run is the mass of the produced hadron, and p q T is the 
transverse momentum of the light quark. From above, we 
get the form of light quark fragmentation function: 

D f {z, E q ) = a(E q - E )(l - m 2 H /p 2 qT )P(z, a) , (2) 

where Eq, a and a are parameters to describe the size and 
the shape of the LQF effect. 



3 Calculations for charm hadroproduction 

For high energy hard scattering of hadrons A and B, the 
inclusive cross section for the production of hadron C can 
be factorized as 



E c d 3 a A B- 



■cx 



d*Pc 



E 

abed 



dx l dx 2 f a A {x 1 ,Q 2 )f h B {x 2 ,Q 2 ) 



d&„ 



-D c c {z,Q 2 )—m 



dt 

where fjAxi, Q 2 ) and fg(x2, Q 2 ) are the parton distribu- 
tion functions, with X4 being the momentum fraction car- 
ried by the parton in the infinite momentum frame (IMF). 

is the subprocess cross section, with t being the 



da 



dt 



parton level Mandelstam variable. D^(z,Q 2 ) is the frag- 
mentation function of parton c into hadron C, with z being 
the energy fraction of parton c carried by the hadron C, 
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z = Ec/E c , in the parton cm. frame. Q 2 is the factoriza- 
tion scale, which can be taken as the squared transverse 
momentum of the subprocess. 

Now we consider the process ir~N — > D ± X. Since 
charm hadrons is always regarded as the fragmentation 
product of charmed quarks, with the LO subprocesses 
gg — > cc and qq — > cc, and charge symmetry assump- 
tion £>f + (z, Q 2 ) = DP~ (z, Q 2 ), the produced D+ and D~ 
should be symmetric. However, if the light quark fragmen- 
tation into charmed hadrons is non-negligible, additional 
contribution from the LQF effect should be considered ac- 
cording to Eq. © • The largest contributions from the LQF 
effect_are d + g — ► d + g, with d —>_D~ , and d + g — > d + g, 
with d — > D + . Since more d than d quarks exist in 7r~ and 
the nucleoli, more d than d quark will contribute to the 
process, and more D~ than D + will be produced, just as 
experimentally observed. 

LO PQCD calculation of charm hadroproduction for 
fixed target experiments with factorization formula is a 
debated area, since Q 2 is generally only a few GeV 2 in 
these energy region, and higher order corrections could 
be large. The NLO calculation [7] produces an increasing 
factor of about 2 relative to LO calculation, with a simi- 
lar shape for differential cross sections. However, since we 
only aim at a sketchy estimate of the LQF effect in this 
work, and since the asymmetry of Eq. JT]) that we calcu- 
late is a ratio of the cross sections, whose uncertainties 
may be cancelled to some extent, we expect reasonable 
results within uncertainties from LO calculation. 

The differential cross section for the inclusive produc- 
tion of hadron C as a function of xf can be expressed 
as 



da a 



B^CX 



dxf 



J2 J dx 1 dx 2 dzf%{x l ,Q 2 )f b B {x2,Q' 2 

abed 



x%^ c C (z,Q 2 ) 
dt 



dxi 



dt 



(4) 



where xf is the Feynman variable for hadron C, xp = 
2P* /V~S, with P* being the momentum along the incident 
direction in the cm. frame of the interacting hadrons A 
and B, and 5* is the squared center of mass energy S = 
(Pa + Pb) 2 - 

When the masses of quarks and hadrons are neglected, 
Pc = zp c will be hold in any reference frame, and xp can 
be expressed as 



Xp 



(xi + X 2 )t + X\S 



(5) 



different for charmed quark fragmentation and light quark 
fragmentation into charmed hadrons. 

In the case of charmed quark fragmentation, with con- 
sideration of charmed quark mass m c and the produced 
charmed hadron mass to#, Eq. (J5]) and Eq. © can be 
corrected as 



xf = z 



Xi-X 2 ,X X +X 2 ,~ 2 , s x ,-s, 

7, + z : (*-™*c + o)' ( 7 ) 

Is I 



and 



dxi 



dt 



X1X2S 
z'(x\ + x 2 ) 

e c ), with eu = 



(8) 

2m 2 H /{z 2 s) and 



where z' — z(l — en 
e c = 2m 2 c /s. 

In the case of light quark fragmentation, when the pro- 
duced charmed hadron mass mu is considered, and notice 
that the fragmentation function is defined in the nucleon 
rest frame with z being the energy fraction of the produced 
hadron, Eq. ([5]) and Eq. ([6]) can be corrected as 



xf 



{xi + x 2 )i + XiS 



2ze^ EA 



M)p q 



t H- 



and 



dxF 


-1 


di 





x\ + x 2 ze^xiEA ,2p, 



Ms 



(9) 



,(10) 



where Ea is the incident energy, which is 500 GeV for the 
7r~ beam in the E791 experiment, and M is the mass of the 
nucleon as target; = m 2 H / '(2z 2 E 2 ) , and E q and p qz are 
the energy and longitudinal momentum of the outgoing 
light quark in the nucleon rest frame, 



E q = x x E A {l 



t , 



x 2 M 



. , t x 2 M 
Pqz = [XiE A + x 2 M) — + x x E A H — 



(11) 



(12) 



Similarly, the differential cross section as a function of 
can be expressed as 



do ab- 



■cx 



dP 2 



= ^ J dxidx 2 dzf A (xi,Q 2 )fg(x 2 ,Q 2 ) 

abed 



da a 



dt 



D^(z,Q 2 



dt 



.(13) 



where s — x\x 2 S is the squared center of mass energy of 
the subprocess. Thus the term \ in Eq. (|4]) is 



When the mass effect is neglected, P% and 
expressed as 



can be 



dxi 



dt 



x\x 2 S 
z{x\ + x 2 ) 



(0) 



However, the masses of c quark and D ± meson should 
be considered in our case, and their effect can be taken 
as corrections to Eq. and Eq. ((6|). The corrections are 



and 



p2 



■«7 + t > 



dt 



z 2 \s + 2t\ 



(14) 



(15) 
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In case of charm quark fragmentation, with the inclusion 
of the mass effect, Eq. ([T4| and Eq. (fT5|) can be corrected 
as 



and 



4P| 



dt 



z' 2 \s + 2(t- 



mt 



(16) 



(17) 



where z' = z(l — en + e c ), with en = 2m 2 H /(z 2 s) and 
e c = 2m 2 /s. In case of light quark fragmentation from 
the nucleon rest frame, with the mass of the produced 
hadron m# being considered, Eq. (|14p and Eq. (fTS"|) can 
be corrected as 



,t z 



and 



P 2 = -z 2 (l-2e-)(j+t) 



dP 2 



(18) 



dt 



2e%) + 2t\ 



where = m 2 H / '(2z 2 E 2 ) , and i£ 9 is expressed in Eq. (fTTj) . 

In our calculations, s > (2mn) 2 and z > 2mjj/vi are 
required in cc production and fragmentation processes, 
and mjj /E q < z < 1 — {m,A a — M) / E q is required in light 
quark fragmentation. To compare our calculation with the 
E791 data, restriction —0.2 < xp < 0.8, just as in the ex- 
periment, is obliged in the calculation of P^ distribution. 
As to the calculation of xp distribution, P^ is indirectly 
restricted by Q 2 > Qq. 

We use CTEQ6L1 parton distributions [29] for the nu- 
cleon, and we employ the LO parametrization forms of 
Ref. [30 for the parton distributions of the 7r~. Since 
the E791 experiment uses a target mainly of carbon, we 
take it as an isoscalar target for the nucleon parton dis- 
tributions. We consider LO subprocesses gg — ► cc and 
qq — > cc for charm quark production, and various LO 
subprocesses including qg — > qg, qq' — > qq' , qq — > qq^ 
qq — ► qq, qq — ► q'g', — > qg, for light quark d and d 
production. The cross sections ^ for these subprocesses 
can be found in Refs. [3T1I321I33] . We use the LO running 
a s with Aqcd = 215 MeV for n/ = 4 and vIqcd = 165 
MeV for nj — 5, as specified in CTEQ6L1. We set factor- 
ization scale Q 2 = p 2 T + m 2 for charm quark production 
processes, and Q 2 = p 2 T for light quark production pro- 
cesses, with m c = 1.5 GeV and Q 2 , = m 2 H w 3.5 GeV 2 , 
where p,jT is the transverse momentum of quark q. We 
use Peterson parametrization [34] for charm fragmenta- 
tion function Dj?(z, ep), with ep — 0.08 and the fragmen- 
tation fraction 0.23 for c — > D + (including D + from D* + 
decay) [35||36]. For the light quark fragmentation d — > 
and e? — > _D + , we use the parametrization of Eq. @ and 
tuning the parameters to describe the observed asymme- 
try. 



4 Results and discussions 

From numerical calculations, we can find good descrip- 
tions of the observed asymmetry from the E791 experi- 
ment. Fig. Q] shows a set of the best fits for both xp and 
P^ distributions with the following parameters, Eq = 40 
GeV, a = 2.0 x 10~ 5 GeV -1 , and a = 0.6. Notice that our 
A(xp) result is consistent with the data of all xp region, 
while A(Pj-) result could not give a good description of 
the data below P^ ~ 2 GeV 2 . This discrepancy, however, 
can be attributed to the unaccounted small random trans- 
verse momentum of the produced hadron relative to the 
direction of the light quark when it fragments. When this 
effect is considered, the sharp peak of D produced from 
LQF at low transverse momentum could become broader, 
and the asymmetry curve should get greatly balanced at 
low PjL Meanwhile, the general features of A(xp) and 
A(P%) (P| > 2) do not change. 

Eq is a parameter that describes an energy threshold 
for light quark fragmentation into charmed hadrons. There 
are still uncertainties in its value, as our result is insensi- 
tive to this. Nevertheless, the overall size of the LQF effect 
(about a(E q — Eq)) is quite stable in our fits, when chang- 
(19) ing E Q from 10 to 80 GeV, to be about 2.4 x 10~ 3 near 
E q = 160 GeV, which is a typical energy under study. 



a describes the shape of the fragmentation function for 
the LQF effect, and a = 0.6 corresponds to a broad peak 
around z — 0.4. 

In this work, we attribute the leading particle effect for 
charm hadroproduction to the LQF effect from the pick- 
up process of the produced light quarks. Nevertheless, we 
do not exclude other possible contributions that lead to 
the fragmentation of light quark into charmed hadrons, 
which should still be explored by some dedicated work. 

From our LO calculation, the LQF effect could give 
a good description of the observed Z) ± asymmetry. On 
the other hand, we can get an estimate of the shape and 
the size of the LQF effect from the fitting results. When 
z is integrated out, the size of LQF effect can be drawn 
from Eq. © to be about a(E q — E ) when p 2 T » m 2 H , 

with E = 40 GeV and a = 2.0 x 10~ 5 GeV" 1 , one gets 
Df(E q ) = 2(E q -40) x 10~ 5 . For a typical value E q = 160 
GeV, one gets Df — 2.4 x 10~ 3 , which is very stable in 
our fit as we have mentioned above. However, there are 
still uncertainties for this fitting result both from LO cal- 
culations and from the parametrization of the light quark 
fragmentation function. Since the energy of fixed target 
experiments is limited, higher order effects will contribute 
and the calculated LO cross sections show a scale depen- 
dence. Meanwhile the low transverse momentum suppres- 
sion factor (1 — m 2 H /p^ T ) in our parametrization gives an 
extra uncertainty by a factor of 1/3. As a sketchy esti- 
mate, we set an increasing factor of 2 and a decreasing 
factor of 3 as its uncertainties around the central value of 
our fit. Thus, for E q — 160 GeV, we have an estimate of 
0.8 x 10^3 <D f < 4.8 x 10~ 3 i.e.Df = (2.4+ 2 ;*) x 10~ 3 
for the size of the LQF effect. 

Now we compare our result with the previous estimate 
of the LQF effect from neutrino induced dimuon produc- 
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Fig. 1. Comparison of the asymmetry variable of Eq. |TJ from the LQF effect of this work with the E791 data [3]. The solid 
lines are our calculation results with parameters Eq = 40 GeV, a = 2.0 x 10 GeV - , and a = 0.6. 



tion process |27j . In the previous work, the light quark 
fragmentation rate, defined as D q = D® + , is differ- 
ent from the definition in this work, which is the rate of 

D f ~ Z>P ' 



d —> D including D from D* 



B*D°' , where and D®* 



decay, i.e., u$ — 
indicate direct fragmentation 
into D or D* meson from light quarks, and B* ^ 1/3 is the 
fraction for D*~ decay to D~ [37] . With an estimate for 
the ratio /D® = 3/1 from spin counting, we find that 
the D q in previous work should be one half smaller when 

D q can be drawn from 



it is compared to Df in this work 
Eq. (13) of Ref. E7], D q /J c » 0.20 /tr^ 



with 



an estimate for f c = 0.86 and taking cr„- „- /cr„- „+ = 
(3.5 ± 1.6)% [T5], which is the prompt dimuon rate for 
100 < E vis < 200 GeV with cut > 6 GeV. From above, 
one can get D q = (6.0 ± 2.7) x 10 -3 . Now we can com- 
pare \D q = (3.0 ±1.4) x lO -3 with D f (Eq = 160 GeV) = 
(2.4^^'g) x 10 -3 , for they have similar energy range. From 
above, we find that the size of the LQF effect from the 
estimate of this work of charm hadroproduction is a lit- 
tle smaller but consistent with that of the estimate from 
neutrino induced prompt fj,~ pT production process. Esti- 
mation of the LQF effect from data in the previous 
work yields a much larger size, which is not supported by 
this work. 

In our previous work [27] , we have shown that the LQF 
effect could influence the measurement of nucleon strange 
asymmetry from charged current charm production pro- 
cesses. When the LQF effect is large enough, it can balance 
the effect from nucleon strange asymmetry, which could 
explain why CCFR and NuTeV experiments do not show 
evidence for the nucleon strange asymmetry. From the es- 
timate of the LQF effect in this work, the size of the LQF 



effect at E a = 160 GeV is about 



io- 



which 



is still twice or more smaller to compensate the predicted 
nucleon strange asymmetry that can explain the NuTeV 
anomaly [21^l23ll24l . However, in CCFR and NuTeV 



experiments (neutrino energy up to 600 GeV), the energy 
of E q can be much larger and the LQF effect is possible 
to be greatly enhanced. Thus in CCFR and NuTeV, the 
LQF effect is still possible to influence the measurement 
of nucleon strange asymmetry to a large extent. Further 
study of the LQF effect under various processes will be 
helpful to further clarify this effect. 



5 Summery 



In this work, we have attempted to explain the leading 
particle effect in charm hadroproduction by the picture of 
light quark fragmentation into charmed hadrons (LQF). 
We can obtain good descriptions of the observed D ± asym- 
metry from the E791 experiment for both xf and dis- 
tributions. Although the uncertainty from our LO PQCD 
calculation is large, we find a total size of the order of 10 -3 
for the LQF effect at a typical energy E q — 160 GeV of 
the light quark. This result is consistent with the estimate 
of the LQF effect from neutrino induced prompt 
data in our previous work. The influence of the LQF ef- 
fect to the measurement of nucleon strange asymmetry in 
neutrino induced charged current charm production pro- 
cesses could be non- negligible. 
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